Fields radiated by a nanoemitter in a graphene sheet 
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Abstract 

The extraordinary properties of graphene make it a very promising material for use in optoelec- 
tronics. However, this is still a nascent field, where some basic properties of the electromagnetic 
field in graphene must be explored. Here we report on the fields radiated by a nanoemitter lying 
on a graphene sheet. Our results show that this field presents a rich dependence on both frequency, 
distance to the source and dipole orientation. This behavior is attributed to distinct peculiarities 
on the density of electromagnetic states in the graphene sheet and the interaction between them. 
The field is mainly composed of an core region of high-intensity electromagnetic field, dominated 
by surface plasmons, and an outer region where the field is practically the same it would be for 
an emitter in vacuum. Within the core region, the intensity of the electric field is several orders 
of magnitude larger than what it would be in vacuum. Importantly, the size of this core region 
can be controlled thorough external gates, which opens up many interesting applications in, for 
instance, surface optics and spectroscopy. Additionally, the large coupling between nanoemitters 
and surface plasmons makes graphene sheets a propitious stage for quantum-optics, in which the 
interaction between quantum objects could be externally tailored at will. 

PACS numbers: 42.25.Bs, 41.20.Jb, 42.79.Ag, 78.66.Bz 
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Introduction. - Graphene has attracted recently a great deal of attention due to its amazing 
electronic properties [TJ On top of fundamental issues, such as being the thinnest pos- 
sible two-dimensional electron gas (2DEG) and possessing charged massless quasiparticles, 
graphene presents very interesting material properties for its use in electronics. Specially 
noteworthy are both the possibility of controlling the chemical potential, /J (and thus the 
conductivity a) through gate voltages and its high mobility at room temperature. Also the 
optical properties of graphene are exceptional [H 0], as they allow the visualization of a 
material that is just one atom-thick, and even differentiate regions with different number of 
atomic layers [3 E] • Remarkably, even one-atom thick graphene can bind electromagnetic 
(EM) modes, despite being almost transparent. These surface EM modes are surface plas- 
mon polaritons (SPP) with a transverse-magnetic (TM) polarization [7HTU] when Imcx > 0, 
which occurs below a critical frequency uq that depends on /x, thus being externally tunable. 
For frequencies above u , Imcx < and graphene supports transverse electric (TE) bound 
EM modes [TT], which are reminiscent of the guided modes in a dielectric film. The existence 
of SPPs in graphene, with their similarities to those supported by metal surfaces, brings to 
the fore the possibility of using graphene for many of the functionalities sought within the 
field of Plasmonics. For instance, graphene could be used to built terahertz (THz) switches 
[12], and a graphene sheet with a spatially non-uniform conductivity is a promising platform 
for both THz flatland metamaterials and transformation optics [13J. 

For all these optical functionalities, an essential requirement is the efficient excitation 
of the surface EM modes. In this work we analyze how these bound modes are excited 
by an emitter of deep subwavelength dimensions (quantum dots, molecules, dielectric or 
metallic protuberances, etc.) in graphene. Additionally, we study the EM field pattern 
at the graphene layer, which is an essential ingredient for the understanding of the EM 
properties of nanostructures placed on a graphene sheet and the effective EM interaction 
between them. Our results show that the field patterns are drastically different for different 
frequency ranges and orientation of the dipole. While a nanoemitter hardly couples to bound 
TE modes, SPPs in graphene can be efficiently excited in a frequency range that depends 
on the chemical potential of the 2DEG, which allows the control of their excitation by 
external gates. This suggests graphene as an excellent candidate for supporting adressable 
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FIG. 1: (a) The geometry of the studied system, (b) Dispersion relation q(p,) (continuous curves) 
and propagation length (discontinuous curves) of SPPs and TE bound modes, where q is the 
wavevector in units of uj/c and Vt = hoj/fi. These magnitudes are obtained for the case in which 
fx = 0.2eV (48 THz) and T = 300K. The inset shows a = 2ira/c in units of the fine structure 
constant ao ~ 1/137. 

long-distance entanglement between qubits |14j . 

Description of the system under study.- We consider a point emitter (placed at the ori- 
gin, with dipole moment p) lying on a free-standing graphene sheet (covering the plane 
z = 0). We will concentrate on the study of the electric field E at the graphene sheet, 
which is the relevant quantity concerning the coupling of subwavelength objects placed on 
graphene. Graphene is represented by its in-plane complex conductivity a, that is a function 
of both frequency, u, and material properties, such as chemical potential /x and tempera- 
ture T. In this work, we consider T = 300K and /i = 0.2eV (48THz), which are typical 
experimental values. The conductivity a is taken from calculations based on the random- 
phase-approximation [ToT[T7] . and is represented in the inset to Fig. lb. For this set of 
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parameters the imaginary part of the conductivity vanishes at Hooq = 0.326^. Nevertheless, 
we have checked that all results presented in this paper are qualitatively valid for other 
values of \x (provided that u is scaled in the same manner) and T. We have also checked 
that our results hold true under the presence of a dielectric substrate, if the substrate does 
not support optical resonances in the frequency range of interest. The interaction of EM 
modes in graphene with those of a polar substrate is beyond the scope of this paper and is 
left for future work. 
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FIG. 2: The different components of the angular spectrum Dyadic D as a function of the normalized 
in-plane wavevector q = ck/u. The main figure corresponds to the graphene sheet (for the case 
T = 300K and fx = 0.2eV) while the inset is for vacuum. In both cases f2 = 0.02, which corresponds 
to v ~ ITHz (A = 0.3mm) for the [i considered. 

Mathematical formulation.- The electric field radiated by the dipole can be written as 
E(r) = G(r) p, where G(r) is the Green's dyadic of the problem. We express distances R 
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in dimensionless units as r = k u R = 2tiR/\, being k u — co/c, c the speed of light and A the 
wavelength of the EM field. Similarly, we define dimensionless quantities for wavevectors 
(through q = k/k u ), conductivity (a = 2na/c) and frequency (ft = ftw/n). Due to the 
symmetry of the problem, it is convenient to work in cylindrical coordinates, r = (r, 9) (see 
Fig. 1). Following standard techniques in electromagnetism (see for instance [21]). G can 
be represented in terms of Sommerfeld integrals: 

■ik « r°° 

G ^ = ^ E / d< i D lM) J iM r ) (!) 

t=TE,TM ^° 

where i,j are r, 9 or z, J[j(qr) are combinations of Bessel functions and D\- are the compo- 
nents of the angular spectrum dyadic (see the Supplementary Material for the derivation and 
explicit expressions). The non-zero components of DJj are represented in Fig. 2 for both the 
graphene sheet (main figure, for Q = 0.02) and free space (inset). For this frequency, SPPs 
are supported by the sheet, as clearly seen in the presence of a pole in the TM components 
of the angular spectrum of graphene at q p = k v jk^ ~ 1.9. The location of poles in complex 
q— space, q p , provide the spectral region of existence of SPPs and bounded TE modes, and 
their corresponding wavevector, Re(g p ), and propagation length, L p = A/(27iTm(g p )), which 
are rendered in Fig. lb. 

In principle, Gij(r) must be computed numerically [20]. However, some heuristic proper- 



ties can be extracted from the form of Eq, A. 16 and with the visual help of Fig. 2. Equation 



1 shows that the EM field originates from the interference of a continuum of oscillatory 
functions, JJj(qr) (which are given by the symmetry of the problem), each of them launched 
with an amplitude DJj (which depends on the material properties). For large distances, 
the oscillatory functions oscillate very quickly with q, and the contributions to the integral 
of "smooth" (in a scale ~ 1/r) parts of the angular spectrum cancel out. Thus, only the 
sharpest spectral features in DJj are relevant at long distances from the source. These can be 
divided in three classes: (i) square-root singularities, which appear for some components at 
the light-cone (q = 1) and lead to a 1/r dependence of the field at large distances, (ii) poles 
associated to bound EM modes which decay exponentially with distance both in the normal 
and in-plane directions and (iii) kinks (points where the derivative of the angular spectrum 
is discontinuous) which lead to 1/r 2 decays. The expressions for these contributions are 
given in the Supplementary Material. As we will show later, all these spectral features show 
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up in the EM field in different spatial regions. 

In what follows we present results for the electric field at the graphene sheet radiated 
by a nanoemitter as a function of frequency ou. Let us analyze separately the emission by 
considering the different components of the dipole. 

Perpendicular dipole, p = pe z .- Figure 3 presents the spatial dependence of the electric 
field at the graphene sheet. Due to the symmetry of the problem the field does not depend 
on the polar angle 9, as illustrated in panel (a). Panels (b)-(d) render the results for four 
different frequencies, each one being representative of a different regime. In each case, the 
exact electric field (E, red lines) is plotted together with the SPP contribution (E SPP , black 
lines) and, for reference, the electric field that the dipole would radiate in free space (FS) 
(E FS , blue lines). 

Panel (b) renders the results for = 0.02, that corresponds to a frequency much smaller 
than the chemical potential. In this case, k p ~ k u and L p ^> A. As the figure shows, the 
exact field virtually coincides with the SPP contribution even at distances much smaller 
than the wavelength. At larger distances the exponentially-decaying SPP dominates until 
it is superseded by a algebraically-decaying field, which turns out to be well approximated 
by the free-space radiation (see inset to Fig 2b). The crossover distance R c typically occurs 
at R c ~ 10Lp which, for this frequency regime, may be hundreds or even thousands of 
wavelengths. The field oscillates with a similar period before and after R c , presenting an 
interference pattern at distances close to R c , where SPP and FS fields are of the same order. 

Figure 3c is representative of an intermediate frequency range (Q = 0.15), in which 
k p ^> k u and L p is a few times the wavelength. In this range the field can also be well 
approximated by the sum of two terms: the SPP (which dominates at short distances) and 
a "free-space" wave (prevailing at long distances). However now the spatial period of the 
EM field at R < R c and R > R c are vastly different (see insets to Fig. 3c). What is 
important to note is that now the E-field amplitude in the close proximity of the emitter is 
much larger (around one thousand times) than that obtained for free-space radiation. For 
higher frequencies, L p < A and the SPP is an overdamped wave (0.4 < il < 1.6 for the 
temperature and chemical potential considered in this paper, see Fig. 1). Still the field at 
the graphene sheet can be represented by a SPP at R < R c and as free-space radiation for 
larger distances (see Fig. 2d for f2 = 1). However, in this case, R c is in the subwavelength 
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FIG. 3: Spatial dependencies of the electric field radiated by a dipole perpendicular to the graphene 
sheet: p = pe z . In (a) the field colormap for log(|E|) is shown for £7 = 0.02. The white discon- 
tinuous line represents the section along which all the ID field dependencies are shown. In (b), 
(c), (d) and their insets, both log(|E|) and Ke(E z ) are plotted as function of distances for different 
Q. In all the panels together with the full field, both surface plasmon and free-space contributions 
are shown. All fields are normalized to k w p/2, which is the largest value of the field radiated by a 
dipole moment in free space at distance r = A. 

regime, and in fact it vanishes as Q — > Qq (Qq = 1.6 for the material parameters considered 
in this paper). Nevertheless, note that the huge amplification of the E-field in the vicinity 
of the emitter is still maintained thanks to the SPP excitation. For frequencies Q > Qq, the 
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graphene layer does not support SPPs but it does support TE surface waves. However, as 
represented in Fig. 2d for Q = 2, the fields at the graphene sheet are very approximately 
the same as if the graphene layer were not present. The TE field is so weakly bounded that 
it virtually does not couple to the point emitter. 

An important message that can be extracted from Fig. 3 is that in all cases where SPPs 
are supported by the graphene sheet, the intensity of the field at distances from the source 
smaller than R c is orders or magnitude larger than in free-space. Our results show that this 
enhancement can be as large as 10 5 — 10 6 . This amplification, combined with the tunability 
of the chemical potential (and thus of O) through doping levels and/or external gates may 
have important practical applications. For instance, whether a subwavelength region close 
to the emitter (and the extent of this region) is strongly illuminated or not can be externally 
controlled by DC voltages. For this, the chemical potential of the graphene sheet must be 
changed so that for the resulting Q = fvjj/fi SPPs are either strongly overdamped or not 
supported at all (see Fig. 2d). 

Parallel dipole, p = pe x .-. Let us consider now a dipole parallel to the graphene sheet, in 
a direction that we take at the x— direction (i.e. 6 = 0). The electric field at the graphene 
sheet is E(r, 9) = (G rr (r) cos 9, Ggeir) sin 9, G zr {r) cos9) T p, where T stands for transposition. 

In contrast to the case of the perpendicular dipole, the electric field radiated by the 
parallel dipole is strongly angular-dependent, and this dependency drastically changes with 
Q, see colormaps in Fig. 4. Overall, the field pattern follows the rules previously discussed 
for the vertical dipole: (i) For Q < Qq, SPPs are sustained and dominate the field in an 
inner core region, with a size related to the SPP decay length and, thus, strongly frequency 
dependent. Outside this core region, the field at the graphene sheet is very approximately 
equal to the one in free-space (ii) For Q > Qq, the core region has vanished and the field 
is practically the same as for free-space. However, in the case of the in-plane dipole there 
are two main differences. The first one is that, strictly speaking, the long-distance limit of 
the in-plane field is not the transversal free-space field. Instead, the asymptotic field arises 
from the kink in the angular spectrum (see Fig. 2), being reminiscent of the Norton wave 
found in the studies of radio transmission at the earth surface [22] and radiation by holes 
in metals [15] . Our calculations show that the Norton wave in graphene is predominantly 
TM polarized. Figure 5 renders a "phase-diagram" showing the main character of the field 
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FIG. 4: The spatial dependencies of the electric field radiated by a dipole parallel to the graphene 
sheet p = pe x for different Q: U = 0.02 (y = ITHz, A = 0.3mm), Q = 0.4 (u = 19.4THz, 
A = 15.5/xm), U = 2 [y = 97THz, A = 3.1/mi). In (a), (c), (e) the field colormaps for log(|E|) 
are shown. The white discontinuous lines represent the directions along which all the ID field 
dependencies are shown in (b), (d), (f). In all the panels together with the full field, both surface 
plasmon and free-space contributions are shown. The normalization of the fields is the same as in 
Fig. 3. 

as a function of both frequency and distance to the source for the two directions of the 
point dipole considered in this paper. Notice however that the small value of the field in the 
region where NW dominate most probably prevents its experimental detection. The second 



9 



difference is related to the polarization of the inner core and outer region. As SPPs have a 
longitudinal field component (thus vanishing at 9 = 7r/2, see Fig. 4a) and FS radiation is 
transversal (vanishing at 9 = 0, see Fig 4e), the polarization patterns of the inner core and 
the outer region are rotated by 7r/2, as illustrated in Fig. 4c. Correspondingly, for Q < Qq, 
the inner core region extends further away at 9 = than at 9 = 7r/2 (see Fig. 4d). 




distance/wavelength 

FIG. 5: Phase-diagram showing the main character of the field at the graphene sheet as a function 
of source frequency and distance to the source. Panel (a) is for a dipole perpendicular to the sheet, 
while panel (b) is for a parallel one. In both cases T=300K and [i = 0.2eV. 

Note added in proof: during the process of submission of this manuscript we have become 
aware of the related work F. H. L. Koppens et al, arXiv:1104-2068vl (2011), where the 
Purcell factor of the point emitter close to graphene structures has been studied. 
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Supplementary Material: Explicit expressions for the conductivity and the Green's 
Dyadic of the graphene sheet 



The model for the conductivity of graphene 



The conductivity of graphene has been computed within the random phase approximation 
[ToT - HT] . In terms of the chemical potential \i and the temperature T the conductivity can 
be expressed as: 

G 'intra Winter ("^-"^) 

where, for relaxation times much larger that w -1 , the intraband and interband contributions 
are: 



2ieH 

@intra ~ ~ m 



2 cosh 



Winter 



e 

Ah 



2 r 



1 1 [Vt-2 
— I — arctan 

2 7T \ 2t 



In 



(A.3) 



2tt (fi - 2) 2 + (2t) 2 _ 

In this expressions Q = hw/fi and t = T/fi, with T expressed in units of energy. 

A finite relaxation time, r, can customarily taken into account by substituting u> by 



UJ + IT in a intra . 



The electromagnetic Green's Dyadic for a 2D electron gas 




FIG. 6: The geometry of the problem: the dipole with the dipole moment p placed at the origin 
creates an electric field E at the point (r, 6, 0). The plane of 2DEG coincides with z = 0. Both 
cartesian and cylindrical coordinate systems are shown. 

Throughout this Supplementary Material, we express coordinates in the in-plane (R) and 
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normal (Z) directions to the graphene sheet in dimensionless units as r = k^H = 27rR/A 
and z = k^Z. 

The electric field E(r, z) emitted by an electric dipole, with dipole moment p(z') and 
placed at the point (r' = 0,z'), is given through the Green's dyadic (GD) G(r, z, z') = 
G(r, z\ r' = 0, z') by the following relation 

E(r,z) = G(r,z;z')p{z'). (A.4) 



General form of the Green's Dyadic under the presence of a 2D electron gas 

The GD can be expressed (see for instance [191 12] ) i n terms of the electromagnetic 
modes in free-space u qr e* qr+i,?z2 , characterized by their in-plane momentum q = \s./k u and 
polarization r = TE, TM): 



G(r, z\ z') = G (r, z; z') + G R (r, z\ z') 



(A.5) 



with Go being the GD in free space, 

G (r,z;z') = J2j ^<r<^^\ 
and Gr being the contribution due to the reflection from our 2D system 



(A.6) 



G R (r, z; z') = J2j 2y r 9 U q- u 



(A.7) 



In these expressions "+" corresponds to the regions where z > z' while "— " is for z < z', 
and the unitary vectors characterizing the polarization of each mode are: 



u 



qTE 



' u qTM 



(A.8) 



In the above expressions q z = y 1 — q 2 is the normalized z-component of the wavevectors 
and r r is the reflection coefficient for the 2D electron gas. These coefficients are obtained so 
that the magnetic H and electric E fields satisfy the following continuity relations: 



e z x (E_ - E+) 
e 7 x (H - 



0. 



Aw. An „ 
—j = — ae z x E H 

c c 



(A.9) 
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where E_ (H_) and E + (H+) stay for the electric (magnetic) fields in the regions of negative 
and positive z, respectively, and e z is the unitary vector along the +z direction. As a result 
of the matching we have 



„TE 



—a 



TM _ ~ a Qz 



q ot + q z ' q aq z + 1 ' 

with a = 2ira/c, being the dimensionless 2D conductivity. 
Explicitly, we have for G = G% M + Gq E 



iqr+iq z \z-z'\ 



87T 2 J q z q 



-QxQy ql 




"ik-w f dq jqr+iq z \z-z'\ 



^7T 2 J q 2 



( 2 2 \ 

QxQz q x q y qz Tq x q x 

q x q y q z q 2 y q z Tq y q 2 

\rq x q 2 Tq y q 2 q 4 /q z J 



Analogously, G R = G T R M + G T R E reads 



R 8vr 2 J q z q 2 q 



^ q 2 y -q x q y 

-q x q y ql 




V 



G™(r) 



ik^ f dq TM - +iqz(z+zl) 
8tt 2 J q 2 q 



( 2 2 \ 

' q x qz q x q y q z q x q x 

q x q y q z q 2 y q z q y q 2 

\-q x q 2 -q y q 2 -q 4 /q Z j 



(A.10) 



(All) 



(A.12) 



Green's Dyadic for the dipole placed at the interface z = 



For the case when the dipole is placed onto the graphene sheet, that is when we take 
z' = + in the previous expressions, the Green's dyadic G(r, z) = G(r, z, z' = + ) can be 
greatly simplified. 

The symmetry of the problem makes it convenient to work in cylindrical coordinates 
(r, 6, z), see Fig. [6 



x = rcos6, y = rsin9, z = z. 



(A.13) 
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In this system of coordinates the Green's dyadic can be obtained from the one in cartesian 
coordinates through 



Qcyl 



indrical 



T G 



1 /^cartesian 



T. 



(A.14) 



T 



V 



A 



(A.15) 



where 

cos 9 — sin 9 
sin 9 cos 9 
II 

In the polar coordinates the Green's Dyadic can be expressed in the form of Sommerfeld 

integrals involving the Angular Spectrum Dyadic D(q): 

i_K 
8tt 



Gij 



r=TE,TM 



(A.16) 



with 
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J_(gr) Ji(qr) 
J" 1 {qr)= ./.U/r) • (A.17) 
Jx{qr) Jo(qr) , 
In this expressions, J±(qr) = Jo(qr) ± ^(qr) and J n (qr) is the Bessel function of nth 
order. The non-zero components of the angular spectrum dyadic are Dj r E = Dj E = t^ E q/q z , 



D 



TM 



D 



TM 



t™nn D 



TM 

rz 



D 



TM 



-2iq 2 t™/q z , and D 



TM 

zz 



2qH T q M /q z 



SPP contribution to the field 



The contribution from the SPP is calculated from the residues of the poles in the TM part 
of the Green's dyadic. For the case of a vertical dipole, this calculation has been performed 
in Ref. [101 In order to find the SPP field, all Bessel functions must be expressed in terms of 
Hankel functions of first (H^) and second (H^) kind. Then the contour integral is deformed 
and closed either in the upper half-plane of complex variable q (for terms involving H^) or 
in the lower half-plane (terms where appears). As the SPP pole is enclosed only when 
closing the integration contour in the half-plane, the final result is 
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k u q p 
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^iq^HY 1 (q p r) 2<$a?H™ (q p r) j 
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where H±\qr) = H^^qr) ± H^\qr) and Hn\qr) being the first-kind Hankel function of 
nth. order and q p = \/\ — a~ 2 with q pz = —or 1 are the normalized k-vector components of 
the SPP. 



Asymptotic of the Green's dyadic at the surface z = 0: the leading order 

For large values of r, such that r ^> 1 the asymptotic expansion of GD can be performed. 



As in the previous subsection, the limits of the integral (A. 16) are extended to the whole 
real q— axis using the Hankel functions, and then the long-distance asymptotic value of these 
functions is taken. The result reads 



G(r,z) 
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(A.19) 



dqy/q-e iqr+iq * z 

\ aq z + l ^ q z (aq z + l) J 

On the plane z — 0, the main contribution is expected to proceed from the branch point 
q z = 0. Then the computational procedure consists in transforming the integral to such a 
form so that the integrand is proportional to l/q z . The element zz already contains this 
factor, so that the transformation for it is not necessary. The other elements should be 
integrated by parts. After the transformation the integrand is expanded in the vicinity of 
q z = and then using 



„iqr 



dq— = irH^\qr), 



(A.20) 



the integration is trivially performed: 
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(A.21) 



Finally, substituting Hq (qr) by its asymptotic value, we arrive at 
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(A.22) 
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The dependency ~ 1/r is exactly the same as for the case of free space (spherical- type 
wave in vacuum) while the ~ 1/r 2 decay is redolent of the Norton wave arising in the radio 
transmission at the earth surface |22j and radiation by holes in metals |18j . 



Analytical form of the Green's Dyadic of free space for z = 



The free space GD follows from Eq. (A. 5) just by setting a = (in which case, the 



reflection coefficient vanishes) so that we simply have G = G . Gq can be rearranged as 

dq 



Go(r,z) 



ik„ 



-q x q y 1 - q* -q y q z 

\-q x q z -q y q z 

Using the Green's function for a scalar potential 



Jqr+iq z z 
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Eq. (A.23) can be written in a differential form 

/ 



G (r,z) 
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(A.23) 



(A.24) 



(A.25) 



After performing the derivatives in (A.25), the result in cylindrical coordinates is: 



^cylindrical^ 



47T r 3 
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(A.26) 
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